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We demonstrate that for general conformal field theories (CFTs), the entanglement for small 
perturbations of the vacuum is organized in a novel holographic way. For spherical entangling regions 
in a constant time slice, perturbations in the entanglement entropy are solutions of a Klein-Gordon 
equation in an auxiliary de Sitter (dS) spacetime. The role of the emergent time-like direction in 
dS is played by the size of the entangling sphere. For CFTs with extra conserved charges, e.g., 
higher-spin charges, we show that each charge gives rise to a separate dynamical scalar field in dS. 


Introduction and summary.- Understanding the struc¬ 
ture of quantum entanglement has emerged as a cen¬ 
tral question in elucidating novel emergent phenomena 
in complex quantum systems. This issue arises in a wide 
variety of research areas, ranging from condensed matter 
physics to quantum gravity. In the former, entanglement 
entropy distinguishes exotic phases of matter, such as 
quantum Hall fluids or spin liquids, while in quantum 
gravity, entanglement plays a key role in the emergence 
of quantum spacetime and also of the gravitational equa¬ 
tions of motion m- A great deal of recent progress in 
this area has come from the AdS/CFT correspondence. 
In this framework, entanglement entropies in the bound¬ 
ary CFT are encoded holographically in terms of the 
Bekenstein-Hawking entropy of extremal surfaces in the 
dual bulk spacetime Ii[7]. 

In this paper, we show that the entanglement structure 
of any CFT has a novel holographic description in terms 
of an auxiliary de Sitter (dS) geometry. As is typical, 
scale in the CFT emerges as the extra holographic direc¬ 
tion, but in the present construction the scale plays the 
role of time in the dS geometry. This structure is rem¬ 
iniscent of the effective causal structure inherent to the 
Multiscale Entanglement Renormalization Ansatz BE 
and the interpretation of these tensor networks within 
AdS 3 /CFT 2 correspondence [TD]. Quite remarkably, we 
find that the emergent dS geometry is the arena for a 
local dynamics, in which one of the degrees of freedom 
is identified with the perturbations of the entanglement 
entropy (EE). We emphasize though that our construc¬ 
tion is not related to the standard AdS/CET correspon¬ 
dence and applies for any CFT in any number of dimen¬ 
sions, without any requirement of a large central charge 
or strong coupling. 

We begin with a d-dimensional CFT in its vacuum 
state in flat spacetime. Now consider evaluating the EE 
for a (d — l)-dimensional ball B of radius R centered at 
X on a fixed time slice. As we review below, for weakly 
excited states, the change in the EE is fixed by the ex¬ 
pectation value of the energy density: 

SS{B) = 27r ' (1) 

While this result is now fairly well known, it went 


completely unnoticed that the integration kernel in the 
above expression is a boundary-to-bulk propagator in d- 
dimensional dS geometry where the radius R plays the 
role of the time-like coordinate 

ds^ = ^ {-dR? + dx^) . (2) 

Hence, 5S{x,R) obeys the Klein-Gordon equation 

(V,V“-m2) dA = 0, (3) 

in this auxiliary dS, where the mass is given by 

rn^L^ = -d. (4) 


This result is the focal point of the Letter. Eurther, as 
we demonstrate below, for CFTs with extra global {e.g., 
higher-spin) charges, there exists one additional dynam¬ 
ical field in dS for each charge. 

With Eq. (|^, the asymptotic boundary data {i.e., the 
behavior at i? = 0) is the expectation value of the energy 
density {Ttt), which sets 5S at very small scales. Then 
the EE perturbations at larger scales are determined by 
the local Lorentzian propagation into the dS geometry. 
Hence the EE for small excitations of the vacuum state 
in any CFT is organized with respect to scale in a novel 
Lorentzian holographic manner. As we discuss below, 
the choice of the asymptotic boundary data implicit in 
Eq. § is precisely that needed to remove the unstable 
modes associated with the mass term 0 being tachyonic. 

We would like to point out that in the AdS/CFT 
framework, the same wave equation (|^ implicitly ap¬ 
pears in mm but their derivation relied on the use of 
holographic EE and the Einstein equations in the bulk. 
However, we reiterate that our construction is not con¬ 
nected to the AdS/CFT correspondence. As we discuss 
below, the present holographic propagation of SS relies 
solely on the so-called ‘first law of entanglement.’ 

Finally, for pure states, the EE of any ball matches that 
of its complement. This condition imposes an antipodal 
symmetry on the solutions of Eq. (|^. Combining this 
property with Eq. Q, we find novel constraints on the 
profile of the energy density — see Eqs. @ and 


First Law as Lorentzian Propagation. - To evaluate EE in 
a quantum field theory, we divide a constant time slice 
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into two parts, a region V and its complement V. Upon 
tracing out the degrees of freedom in U, we are left with 
the reduced density matrix pv describing the remaining 
degrees of freedom in the region V. The EE is then eval¬ 
uated with the standard expression for the von Neumann 
entropy: 

S{V) = -Tr(py logpv) . (5) 


Since the reduced density matrix is both hermitian and 
positive semidefinite, it can be expressed as 


Pv 




tr e“^v ’ 


( 6 ) 


where the hermitian operator Hy is known as the modu¬ 
lar Hamiltonian m- Now for a small perturbation of the 
reduced density matrix py + dp, one can show that the 
change in the entanglement entropy ([^ is given by m 

ss = 6{Hy), (7) 


where S{H) denotes the change in the expectation value 
of the modular Hamiltonian associated with the original 
density matrix py. Eq. 0 is commonly called the first 
law of entanglement, as it is a quantum analog of the first 
law of thermodynamics. 

The cases in which the entanglement Hamiltonian is 
local are rare. One special case is for a GET in its vacuum 
state in d-dimensional Minkowski spacetime and 

where the region of interest is a spherical ball B. In 
this case, the entanglement Hamiltonian takes the simple 
form 

Hb = 27r , ( 8 ) 

where the integral is taken over the ball centered at po¬ 
sition X and with radius R, and Tu is the energy density 
operator. Now, combining this expression with the first 
law Q, we find that for spherical regions, the change in 
the entanglement entropy from the vacuum to weakly ex¬ 
cited states is given by Eq. 0 - Implicitly, we are using 
the fact that the expectation value of Tu vanishes in the 
vacuum. 

As noted above, the integration kernel in Eq. 0 is a 
boundary-to-bulk propagator in d-dimensional dS space 
with the metric 0. Now, one easily verifies that the 
perturbation 5S obeys the wave equation 0 on this aux¬ 
iliary dS. In general, this equation has two independent 
asymptotic solutions which to leading order take the form 

5S ° F{x)/R + f{x)R‘^ + ■■■ . (9) 


Hence, our wave equation admits boundary data with 
conformal weights A = — 1 and d, corresponding to F{x) 
and f[x), respectively. We can then identify Eq. 0 as 
the solution with 

d+1 

F{x) = 0 and /(f) = {Ttt{.x)). (10) 


Time slice I'*' = {a? ] i? = 0} 



FIG. 1. One-to-one mapping between points in dS geometry 
and balls on its future asymptotic boundary I"*". We iden¬ 
tify the latter with the constant time slice in a given CFT. 
The future lightcone of the bulk point intersects on the 
boundary of the corresponding ball. 


Therefore, the expectation value of the energy density 
sets 5S at very small scales {i.e., R 0). Then the 
EE perturbations at larger scales are determined by the 
Lorentzian propagation into the dS geometry, accord¬ 
ing to Eq. 0. Thus, the EE for excited states around 
the vacuum is organized with respect to scale in a novel 
Lorentzian holographic manner. Again, this result ap¬ 
plies for any CFT in any number of dimensions d and 
relies solely on the applicability of the first law. 


We also reiterate that the above choice of boundary 
data (101 precisely removes ‘non-normalizable’ or unsta¬ 
ble modes associated with the tachyonic mass term 0. 

At this point, let us note that the wave equation 0 
is a covariant expression and so it can be rewritten in 
terms of any coordinate system on the dS geometry. 
As usual, changing coordinates in the bulk amounts to 
choosing a new conformal frame in the boundary CFT. 
Hence our holographic construction readily extends to 
the CFT in any conformally flat background. The cylin¬ 
drical background K x is of particular interest be¬ 
low. In this case, a constant time slice corresponds a 
round [d — l)-sphere and the corresponding wave equa¬ 
tion then appears in global coordinates on the dS space, 
e.g., ds^ = L‘^{—dT‘^ + cosh(T)^df2d_i). Explicit exam¬ 
ples of propagation in this cylindrical conformal frame 
and in the flat frame are given in the supplemental ma¬ 
terial, which includes also Refs. da US]. There we also 
give an alternative derivation of the wave equation 0 
based on group-theoretic arguments. 


Auxiliary de Sitter Geometry.- The relation between 
balls on a constant time slice of a d-dimensional CFT 
and the dS geometry is easily inferred as follows: Implic¬ 
itly, we have identified our time slice with the future da 
asymptotic boundary 1+ of dS. Now, for any bulk point 
a: S dS, the intersection of the inside of the future light- 
cone of X with 1+ is then a ball-shaped region, see Fig.[|] 
This establishes a one-to-one map between points in dS 
and balls on the time slice. 
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dB B 1+ = {x B = 0} 



distance’ between these balls as LdR/R, where the con¬ 
stant L becomes the dS radius. This is the only definition 
that respects conformal symmetry on the time slice. 

Finally, we note that the fact that we can associate 
Lorentzian order between balls on a constant time slice 
in Minkowski spacetime does not imply the existence of 
local dynamics respecting this structure. Hence the ap¬ 
pearance of the wave equation ([^ is rather remarkable. 

Antipodal Symmetry.- If we consider excitations which 
are globally pure states, they must satisfy the constraint 
that the EE inside each ball must equal that in its com¬ 
plement. Hence 


SS{B) = SS{B) (11) 


FIG. 2. Penrose diagram of dS. Points correspond to 
whereas horizontal lines represent dB represents the 

spherical entangling surface and B (orange) and B (cyan) 
its interior and exterior. Red lines represent corresponding 
lightcones and x and x their tips, i.e., bulk points in dS corre¬ 
sponding to B and B regions. Note that x is the antipode of x. 


In much of our discussion, the time slice has the topol¬ 
ogy of as assumed in the modular Hamiltonian Q. 

Then this map, and the dS metric in Eq. (i), only cov¬ 
ers half of the dS geometry, i.e., the expanding Poincare 
patch. The missing half can be identified with the com¬ 
plementary exterior regions B on this flat slice. This is 
more easily seen by going to a conformal frame where 
our time slice has topology. As described above, 
this corresponds to choosing global coordinates in the 
bulk dS space. Each spherical entangling surface then 
defines two ball-shaped regions covering complementary 
domains on the These two balls are identihed with 

antipodal points in the dS geometry, see Eig.[^ 

Eurthermore, the Lorentzian structure of the dS ge¬ 
ometry can be given a geometric interpretation directly 
in terms of balls on the GET time slice. Causal relation¬ 
ships between points a; G dS become topological relation¬ 
ships between the corresponding balls Bj,, as illustrated 
in Fig. In particular, we have: 1) yi is in the time¬ 
like future of X <G> By.^ is contained within 2 ) y^ is 
in the null future of x lies within Bx but dBy^ 

is tangent to dB^ at one point; and 3) x and ya (or j/ 4 ) 
are space-like separated <G> the domain of (or By,,) 
extends beyond Bx and vice versa. The latter can be 
refined in a number of ways by taking into account the 
relationship with the antipodal point x, e.g., 3a) x and 2/4 
are space-like separated and 2/4 is in the time-like future 
of X <G> By„ is completely outside of Bx. With d= 2, this 
structure reduces to the relations discussed in Ref. [10] 
for one-dimensional intervals. 

What remains is to give meaning to the scale factor 
of the dS metric in terms of the CFT time slice. This 
can be fixed by considering two concentric balls in 
with radii R and R + dR. We then define the ‘time-like 


In the holographic dS picture, this corresponds to antipo- 
dally even configurations SS{x), i.e., 6S{x) = ^^(x). 

Considering general solutions of the wave equation ^ , 
this antipodal symmetry imposes a relation between the 
two classes of boundary data in Eq. of the form: 


27rV 

rm 




^x' lx — x'P 


fiS!’), (12) 


where the integral is over the future boundary 1+ of dS. 
Now, recall that our entropy configurations 0 are not 
general solutions, but ones with vanishing A = — 1 data. 
Hence combining Eq. ( jlOj ) with the antipodal symme¬ 
try constraint ( | 12 | ), we arrive at a constraint on the en¬ 
ergy density prohle: / d^~^x' |x — x'p {Ttt{x')) = 0. For 
Eq. (Ill to hold for all balls, this constraint must be sat¬ 


isfied for all X. This is equivalent to the vanishing of the 
following moments of {Ttt{x)): 


j d"^ ^x(Tt((x)) = 0, j d'^ ^x X {Ttt{x)) =0, 
and J d'^~^x x^ {Ttt{x)) = 0 . 


(13) 


Note that these constraints can be identified as the van¬ 
ishing of the expectation value of the total energy, the 
boost generators, and the temporal generator of special 
conformal transformations, respectively. Eocusing on the 
vanishing of the total energy, we must recall that we are 
working at leading order in a small perturbation above 
the vacuum. This means that the energy difference will 
appear at higher orders in the expansion. At this point, 
we note that while S{B) = S{B) for all B would certainly 
indicate that the underlying state is globally pure, some 
mixed states will still satisfy the leading order constraint 
0 in our perturbative construction. 

We note that the constraints (13) can also be derived 
directly in the CFT, by comparing Eq. ([^ to the analo¬ 
gous expression for the entropy in the ball’s complement: 


SS{B) = 2-k Id 
Jb 




2R 


(TuiS!)). (14) 
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FIG. 3. The domain of dependence D of a ball B enclosed by 
a given spherical entangling surface. Blue lines are tangent 
to the conformal Killing vector field K^, which gives rise to 
conserved charges (16l. 


Let us mention as well that the constraints (131 can also 


be derived from the results of Ref. [18]. Finally, we in¬ 


clude the analogue of Eq. (131 for the conformal frame 


where the time slice is spherical: 


J (f {Ttt{n)) = 0 , y {Tu{n)) = 0 




Here, points on ^ are parametrized by d-dimensional 
unit space-like vectors n^, with volume form 

Extension to Higher-Spin Charges.- In this section, we 
discuss a generalization of our holographic dS construc¬ 
tion, based on viewing the stress tensor as the special 
spin-2 case of a conserved symmetric traceless current 
d)‘i.../is with arbitrary spin s > 1. The s = 1 case is an 
ordinary charge current J^. The s > 2 case is relevant 
for CFTs with higher-spin symmetry. These include free 
theories in all dimensions, as well as some non-trivial the¬ 
ories in d = 2 {e.g., see [Ullin] and references therein). 


First, let us note that the expression (|^ for the mod¬ 
ular Hamiltonian has a covariant meaning in the CFT 
background spacetime. It is the flux through B of the 
conserved current Jfi ' = T^i,K" where is the (time¬ 
like) conformal Killing vector that preserves the bound¬ 
ary of H, see Fig. This suggests a natural gener¬ 
alization to the spin-s case: for each we de¬ 

fine a charge as the flux through B of the current 
^ . This charge is given by the 

following integral on our chosen time slice 

= (2,)- (^!H!^y"T,,.,(f). 

(16) 

Note that all of these currents are conserved, i.e., 

{s') 

= 0, and hence evaluating this flux through any 


hypersurface which is homologous to B yields the same 
charge The modular Hamiltonian ([^ is the spe¬ 
cial case s = 2, i.e., Hb = One can use the 

new charges to construct new reduced density matrices 
Pb ~ exp [— /isQ*-®^] to refine the measurement of the 

entanglement between B and B. Ref. [2T| studied this 
approach for s = 1, where is an ordinary charge. 
Higher-spin charges QO) with s > 2 were discussed in 
Ref. [221 for two-dimensional CFTs. 


We now observe that the integration kernel in Eq. (16) 


is again a boundary-to-bulk propagator in dS, for a scalar 
with mass given by: 


= —(s — l)(d -I- s — 2). 


(17) 


Thus, the charges can all be interpreted as scalar 
fields in dS, obeying a Lorentzian wave equation with 
mass as in Eq. ( |17[ ). For an ordinary charge with s = 1, 
the bulk field in dS is massless, while charges with s > 2 
correspond to a discrete series of tachyonic masses. Of 
course, these mass values are precisely those that allow 
boundary data with conformal weights d-l-s —2, in agree¬ 
ment with the weights of the densities 


The antipodal symmetry constraints (13) generalize 
most cleanly when recast in their form ( [Ts] ). Then 
the analogous constraint for general spin s is the vanish¬ 
ing of the first s moments of For even (odd) s, 

these constraints lead to antipodally even (odd) wave so¬ 
lutions in dS, so that the charge Q^®^ in every ball equals 
plus (minus) the analogous charge in its complement. For 
example, in the case of s = 1, this reduces to the vanish¬ 
ing of the total charge on the time slice. 

Outlook. - We demonstrated that the conformal group of 
a d-dimensional CFT induces a Lorentzian geometry on 
the space of balls on a constant time slice, which corre¬ 
sponds to dSd. Remarkably, the entanglement of small 
excitations of the CFT vacuum is governed a local wave 
equation (§ on this auxiliary geometry. The expecta¬ 
tion value (Too) is the asymptotic boundary data in the 
dS space, fixing SS at small scales. The EE perturbations 
at larger scales are then determined by the propagation 
into the dS geometry, according to Eq. Hence the 
EE in any CFT is organized with respect to scale in a 
novel Lorentzian holographic manner. We also gave the 
generalization for CFTs with extra global {e.g., higher 
spin) charges, with one dynamical field in dS for each 
charge. 

Our holographic propagation of 6S relies solely on the 
first law of entanglement for spherical regions. The first 
law must apply not just for a particular sphere but for 
all spheres (and their complements) on a given time slice. 
For pure states, the EE of any ball matches that of its 
complement and this introduces an antipodal symmetry 
on the solutions in the dS space. Combining this property 


(111 with the first law, we found novel constraints on the 


profile of the energy density, i.e., Eqs. (13) and (15). 
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Looking forward, it remains to be seen if our holo¬ 
graphic construction can be extended to provide a full 
description of the CFT in terms a local theory of inter¬ 
acting fields, including the metric, propagating in the dS 
spacetime. Such a theory would then provide a novel ex¬ 
ample of the dS/CFT correspondence {e.g., see 
in which the boundary CFT is unitary. 

The present construction is closely related to the pro¬ 
posed description of EE in two-dimensional CFTs in 
terms of integral geometry |10j . Hence, integral geom¬ 
etry may provide an interesting perspective to further 
extend our holographic construction. It may also be that 
our new construction will provide useful insight into ex¬ 
tending the proposal of Ref. m to higher dimensions. 

Of course, a full holographic description would require 
understanding the time dependence of quantities in 
the GET. The natural starting point here would be to 
consider spherical regions not simply on a fixed time slice 
but throughout the d-dimensional spacetime of the CFT. 
The group-theoretic construction presented in the sup¬ 
plemental material suggests that SS now obeys a wave 
equation on the coset SO{2, d)/[SO{l,d— 1) x S'0(1,1)], 
which interestingly has multiple time directions. We are 
currently studying this framework in further detail. 


eCr 


<r„> 



FIG. 4. Top: Rescaled energy density (Tu) for the state (18 ) 
with d — S plotted along the x-axis at t = to-|-5 r, see Eq. (191. 
Bottom: Corresponding change in the entanglement entropy 
SS as a function of dS time R and position along the x-axis. 
The causal nature of propagation is clearly visible. 
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SUPPLEMENTAL MATERIAL 

Explicit Examples of States. - Local propagation of 6S in 
auxiliary de Sitter geometry holds for any states which 
are small perturbations of the vacuum state everywhere 
in space. In particular, we require that for the excita¬ 
tions under consideration, the first law applies not just 
for a particular sphere but for all spheres (and their com¬ 
plements) on a given time slice. Below we discuss two 
examples of such states. 

Let us first consider a pure state in a d-dimensional 
CFT on a plane, which is created by an infinitesimal 


insertion of the energy density operator Tu at time tg+i t 
and position xq 


= (l + eT„)|0), 


(18) 


see also Ref. [19]. The parameter e is taken to be small in 
the sense of e/r'^ ^ 1. The evolution in imaginary time 
T is included to regulate potential UV divergences and 
ensures that the state is a small perturbation of the vac¬ 


uum. The energy density of the state (18) is determined 


by the two-point function of the stress tensor m 


{(l>\Tttit,x)\(l)) = eCr 


(Aa;^ — (At -I- ixY)'^ 


(19) 


(Aa;^-I-(At-I-f r)^)^ l\ 
(Aa:2 - (At-f tT)2)2 ~d)^ 


0[e^ 


where Ct is the central charge and we have defined 
Aa:^ = jaf — xqP and At2 = |t — tgp. An explicit il¬ 
lustration of the energy density for this state and corre¬ 
sponding dS propagation of the perturbation in the EE 
is shown in Eig. Note that the energy density (19) is 
a spherical shell expanding out from (to, Xq) at the speed 
of light. As expected from our general argument, the en¬ 


ergy density profile (19) obeys the constraints (15) and. 


hence, the holographic propagation respects the antipo¬ 
dal symmetry on the auxiliary dS background. 

Our second example is the following mixed state 


p=\0){0\+r,\E){E\, 


( 20 ) 
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where \E) is an energy eigenstate (with constant energy 
density), and ry is a small parameter. In this case, we as¬ 
sume that the constant time slice has topology with 
radius r. Let us now look at (d-2)-dimensional spherical 
entangling surfaces surrounding a cap of the speci¬ 
fied by the angle Oq. The first law reads now 


5S =2tt 


rOo 


r'^-^nd -2 sin ‘^-^0 cte x 


( 21 ) 


cos 0 — cos 00 

X r - X 

sin 00 


77 E 


„d-l 


^d-l 


where f 2 „ = 27r(”+^)'^^/r denotes the volume of a 

unit S". The factors in the integrand then correspond 
to, in order, the volume element of the boundary- 

to-bulk propagator for dS in global coordinates, and the 
(constant) expectation value of the energy density. A 
special case of this expression appears in Ref. [16] , which 
discusses universal thermal corrections to the vacuum 
entanglement entropy. There, the energy is given by 
E = ^, where A is the smallest scaling dimension in the 
spectrum (apart from the identity), and 77 is the product 
of the degeneracy of the energy eigenstate and the corre¬ 
sponding Boltzmann factor, i.e., 77 = Clearly, 


this and other mixed states of the form ( 20 ) violate the 


first constraint in Eq. (15). Hence, the corresponding 6S 


propagates on dS without antipodal symmetry. 

Note that SS in Eq. (21) diverges as Oq —>■ tt, i.e., 


as the dS propagation reaches the past boundary X~. 
This divergence is related to a breakdown of the first 
law and corresponding free propagation in dS space when 
sin 00 ~ rjEr (with 0o > 7’'/2). 

Alternative derivation of wave equation on dS geometry.- 
Let us now present another perspective on the wave 
equation (3). The conformal group relevant for a d- 
dimensional GET is SO{2,d). However, only the sub¬ 
group SO{l,d) leaves a constant time slice invariant. 
Hence the corresponding spherical entangling surfaces are 
mapped onto one another under the action of 50(1,61). 
Now considering the perturbations 6S for these ball¬ 
shaped regions, the 50(1, d) generators Ki act as 


dK,SS [{Ttt)] = -6S [{dK.Ttt)] . (22) 

Here, the on the RHS can be viewed as generat¬ 
ing an “active” conformal transformation that changes 
the GET state, while the dx, on the LHS generates a 
“passive” transformation that instead changes the spher¬ 
ical entangling surface. Comparing now the “active” and 
“passive” action of the quadratic Casimir of 50(1, d), 
= CijdxidKj , we get 

\/^SS [{Ttt)] = SS [(V^Tit)] = -d SS [{Ttt)] , (23) 

where the second expression above uses the linearity 
in Ttt of the modular Hamiltonian ([^. Further, the 
last expression appears after using the fact that the en¬ 
ergy density transforms as a scalar of weight d with 


respect to the 50(1, d) subgroup. Now, a particu¬ 
lar spherical entangling surface is left invariant by the 
stabilizer group 50(1,d — 1). Hence, on the LHS of 
Eq. (23), the nontrivial action of is on the coset space 
50(1, d)/50(l, d — 1). The latter coset is precisely the 
anticipated d-dimensional dS geometry, and becomes 
the d’Alembertian on this space. Hence this group theo¬ 
retic approach produces precisely the Klein-Gordon equa¬ 
tion (3) on the auxiliary dS space. Note that this analysis 
implicitly normalizes the dS radius L to unity. 

Finally, let us mention that the group theoretic argu¬ 
ment above can be also generalized to the higher-spin 
case and, as expected, yields the mass given by Eq. 0 - 
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